Phase diagram of the t-U 2 Hamiltonian of the weak coupling Hubbard model 



00 
O 

o 

(N 



(N 



o 
o 

u 

Q-i! 

03 



o 
o 



> 

On 

m 
o 

oo 
O 



X 



Takashi Yanagisawa° ,h 
a Condensed-Matter Physics Group, Nanoelectronics Research Institute, 
National Institute of Advanced Industrial Science and Technology (AIST) Central 2, 1-1-1 Umezono, Tsukuba 305-8568, Japan 
b CREST, Japan Science and Technology Agency (JST), Kawaguchi, Saitama 332-0012, Japan 

(Dated:) 

We determine the symmetry of Cooper pairs, on the basis of the perturbation theory in terms 
of the Coulomb interaction U , for the two-dimensional Hubbard model on the square lattice. The 
phase diagram is investigated in detail. The Hubbard model for small U is mapped onto an effective 
Hamiltonian with the attractive interaction using the canonical transformation: H e ff = e s He~ s . 
The gap equation of the weak coupling formulation is solved without numerical ambiguity to deter- 
mine the symmetry of Cooper pairs. The superconducting gap crucially depends on the position of 
the van Hove singularity. We show the phase diagram in the plane of the electron filling n e and the 
next nearest- neighbor transfer t' . The d-wave pairing is dominant for the square lattice in a wide 
range of n e and t' . The d-wave pairing is also stable for the square lattice with anisotropic t' . The 
three-band d-p model is also investigated, for which the d-wave pairing is stable in a wide range of 
n e and t pv (the transfer between neighboring oxygen atoms). In the weak coupling analysis, the 
second-neighbor transfer parameter —t' could not be so large so that the optimum doping rate is in 
the range of 0.8 < n e < 0.85. 

PACS numbers: 74.20.-z, 71.10.Fd, 75.40.Mg 



I. INTRODUCTION 

Since the discovery of high-temperature superconduc- 
tors, the strongly correlated electron systems have been 
studied intensively. The effect of the strong correlation 
between electrons is important for many quantum criti- 
cal phenomena such as unconventional superconductivity 
(SC). High-temperature superconductors 0, H, Q as well 
as heavy fermionsjj, d, H, 0| are known as the typical 
correlated electron systems. These systems are modeled 
by the Hamiltonian with the electronic interaction of the 
on-site Coulomb repulsion. Recently the mechanisms of 
superconductivity in high-temperature superconductors 
have been extensively studied u sing the two-dimensional 
Hubbard model [1 ©, 03 El 03 0303 03 03 • 

The superconductivity of the Hubbard model has been 
questioned for many years. It is extremely difficult to 
show the existence of superconducting phase for the Hub- 
bard model in a reasonable way. For the present we can- 
not answer this long-standing question soon. Instead of 
examining the possibility of superconductivity, it is pos- 
sible to investigate possible symmetries of Cooper pairs 
for an effective Hamiltonian with the attractive interac- 
tion. For this purpose effective Hamiltonians have been 
obtained for the Hubbard model. The t-J model is the 
well known effective Hamiltonian derived in the limit of 
the large on-site repulsion U, using the canonical trans- 
formation H t -j = e s He~ s with S oc t/U. On the other 
hand, in the limit of small U, the perturbation theory 
also leads to an effective Hamiltonian with the attractive 
interaction^ El E3 S3, where we have S oc U/t. The 
phase diagram with respect to the Cooper pair symmetry 
can be determined if we solve the gap equation. 

We must notice that we should compare the energy 
with other electronic states to show that the supercon- 



ducting state is indeed stable. For the half-filled band 
with vanishing t' = in two space dimensions, the anti- 
ferromagnetic order parameter for small U is[2l| 
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where c = 3 — It is, however, obvious that the an- 
tiferromagnetically ordered state is unstable away from 
half filling if the Coulomb repulsion U is small. Thus 
we focus on the case of small U for which we have also 
a merit that the gap equation is considerably simplified. 
The purpose of the paper is to determine the gap symme- 
try for the square lattice using the small-t7 gap equation 
derived for the effective Hamiltonian. Although the real 
superconductivity in correlated electron systems should 
be described by a theory of strong-coupling supercon- 
ductivity, the phase diagram can be determined in detail 
using the weak coupling formulation. Precise calcula- 
tions are sometimes not easy at low temperatures in the 
strong-coupling formulation due to the Matsubara fre- 
quency summation and the wave number summation. It 
is important to examine the detailed phase diagram for 
materials belonging to strongly correlated systems such 
as the cuprate high temperature superconductors, the 
organic superconductors, the ruthenate superconductor 
S12RUO4. 

The paper is organized as follows. In Section II the ef- 
fective Hamiltonian is derived using the canonical trans- 
formation. We show that we can derive the attractive ef- 
fective Hamiltonian using some approximations. In Sec- 
tion III the gap equation is shown and the results are 
presented in Section IV. We give a summary in Section 
V. 
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II. EFFECTIVE HAMILTONIAN 

The Hubbard Hamiltonian is 

H = - t E( c ^ +Lc -^ f ' E (c] a c erT +h. c .) 

(ij)a <Cj£>cr 

+ ^ E ™iT™a (2) 



We determine S so as to satisfy H 2 + [S, H ] = 0. We 
find 



S = l y L_ 



— — •4'T c - fc '-9i c -fe-9l c ' £ T- 



Since [S, H 3 ] = 0, we obtain up to the order of U 2 , 



(10) 



where (ij) and <C ji S> denote the nearest-neighbor and 
next-nearest-neighbor pairs, respectively. U is the on-site 
Coulomb repulsion. The unit of energy is given by t in 
this paper. The total number of sites and the number 
of electrons are denoted as N and N e , respectively. The 
half-filled band corresponds to n e — N e /N = 1. 

The effective Hamiltonian is derived using the pertur- 
bation theory for small U . The canonical transformation 
also maps the Hubbard model to an effective Hamiltonian 
with the attractive interaction [2^]. Since no instability 
except superconductivity occurs for small U away from 
half filling, we assume that the pairing interaction is the 
most singular. The procedure of mapping is as follows. 
The Hamiltonian is written as 



where 



H — Hq + Hi + H 2 + H3 , 



ka 



(3) 



(4) 



H ef f = H + H 1 +H 3 + [S, H x ] + - [S, H 2 ] . (11) 



The commutator [S, H 2 ] is evaluated as 



[S,H 2 



E E S kk' 

' kjtk',qjtOpj£p',q'j^O 

x (-5 pfe /c p , 1 .c_ p ,_ 5 ,_ L c_p_,/j i c_ fc ,_ ? _ L c_fc_ g j,Cfcf 
x t t t 

— A t t t 

Op'+q' ,k+qC k i ^C_ k , _ q ^C p ,^Ck\C- p - q i | Cpf 

+ ( Vfc C !'t C -fc'-gX C -fc-<?l C !p'-g' | C -P-<?'1 C PT)' 

(12) 



where 



N 



k^k' 



(5) 



o9 _ 
D kk' — 



- - £,k+ q - ik 



(13) 



Ho 



u 



\ " 



N 2^ c fc'T c -fc'-g| c - fc -?i Cfc T, ( 6 ) 

k^k',q^0 



Since the purpose of this paper is to investigate the pair- 
ing symmetry, we need only the first term and the last 
term of [S, H2]. We find that the average of the second 
and third terms with respect to the BCS wave function 
vanish. Due to the same reason [S, Hi] can be neglected. 
Then the effective Hamiltonian is 



kk' 



The dispersion relation e k for the square lattice is 

efe = —2t(cosk x + cosky) — At' cosk x cosky , (8) 

where /i is the chemical potential. We set = e& — [i. 
Using a canonical transformation, ip = e s ip, we look for 
the solution of the Schrodinger equation H e ffip = Eip. 
The effective Hamiltonian reads 

H eff = eSHe-s = H +[ S,H ]+ 1 - [ S,[S,H ]]+ ... 

= H + Hi+H 2 + H 3 + [S, H + H 1 +H 2 + H 3 ] 
+ hs,[S,H ]] + --- . (9) 



H, 



t-u 2 



ka 



ikc\ a C k a 



u 



Jr 2^ c k>] C -k>l C ~HCk] 
k^k' 



1U*_ 
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E E — — 



^pfe' c P'T C "P-9'i c -p'- 9 'x c -'s-<ji c -fc'-gX c 'sT)' 



(14) 



If we set k = p + q' and p' = k' + q, the first term of 
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[S, if 2] is approximated as 



III. GAP EQUATION 



E E <W 



AT 



p/c' 



x c p'T c -p'-g'l c_ P _< ?'i c -'£'-<;l c_fc_< ?J' C ' £ T 



The gap equation for the t-U 2 model was investigated 
in lief. 1!) . Since the equation was considerably simpli- 
fied for small U, the gap equation was solved without 
numerical ambiguity. We define the order parameter, 



C -k-qj C -k-ql 



Jilo 4 '^ 01 *'"^' ------ --^ 



A fe = Vfefc' (c_ fe /|C fe ' T ). 



(22) 



2 



) E ^'T c -fc'lt fe , 



/- 



k — q 



Using the mean-field theory, the gap equation for the 
Hamiltonian H t _ui is 



x c_ fci c feT 
2 

N 



k^k'-q,q^0 

E 



£k'-q — £,k+q — Zk 



2Ev.i 



(23) 



where Ef. = + A^. We assume the anisotropic order 
parameter given as 



where /fc is the Fermi distribution function, 

1 



(15) 



A k = A • z k , 



(24) 



/fc — 



e /3« fe + 1 ' 



(16) 



where z k denotes the k-dependence of A k . At T = the 
gap equation is written as 



As 



For small [/, the gap equation for anisotropic pairing 
Since the summation is restricted to the small region near is extremely simplified retaining only the logarithmic 
the Fermi surface, we obtain assuming = 



term 1 J 



Ho 



E 



£,k' + k+q - £,(, 



Zk 



= 1 ^(s)° , ^" k+ * : 



Zk' 



(26) 



Similarly the last term of [S, if 2] is written as 

2 



(17) where u>q is the cut-off energy. The critical temperature 
T c is determined by 



H 2b « ( - 



E 



fk' + k+q t t 



-Zk = - 2 , nck'^k' 



2 16 



(27) 



k^-q,k'jt-q 

The resulting effective Hamiltonian is 



(18) for T = T c . For small U, the critical temperature T c is 
extremely small. In this case we can use the following 
approximation, 



H eff = e s He- s = H t _u2 

= E^ fcC ^ Cfc(J +E T/fcfc ' c fc'T c -fe'l c - fe i Cfe T' 



"V^)- tanh(/W2) 



where 

vw = ^ + ^x(k + k'). 

x(fc + fc') is the magnetic susceptibility defined as 

1 fk+k'+q — fq 



(19) 



(20) 



3(w )logwo 



0< ruJ0 







(cosh(/U/2)) 2 



= .g(w )loga;o - 

= 3( w o)l°g^o - 5(0) 

^ f 2e^ 
« .9 log — — 
\7rfc B i c 



/3 c w /2 



cMog ( ^ 







1 



/3 C / (coshx) 5 



(28) 



N q ^k+k'+q 



(21) where we assume that g(t;) is a slowly varying function 
and </(£) is negligible. The equation is written as 



Thus we have reached the effective Hamiltonian up to 
the order of U 2 using the canonical transformation. 
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Since T c is very small, the summation over k' can be 
restricted to the average over near the Fermi surface. If 
we solve the eigenequation 



2 



N 

k' 

the critical temperature is obtained as 

2t 2 

ksTc = 1.13ajQexp 



xU* 



(30) 



(31) 



where the energy unit is given by t. Since A is given as 

2t 2 



A = 2ojQexp 



(32) 



the ratio 2A/fcsT c equals the BCS universal value 
27r/e 7 = 3.53. 



IV. PAIRING SYMMETRY 



A. Method of solving the eigenvalue equation 



< t 



We express zj~ and Vkk' in terms of the polar 
coordinates UM: 



z k = z(£,9), 



(33) 
(34) 



where k is expressed using the polar angle 9: k = (£, 9) 
in terms of the polar coordinates. We consider the gap 
function on the Fermi surface z(9) = z(0, 9). If we define 
x (jS, 9') = x(0, 9, 0, 9'), the gap equation is 

2/ de , p F (6')x{9,6')z{6') = -xz(0), (35) 



where pf(9) is the density of states at the Fermi surface: 

1 , 1 

(27) 



p F (0) = -7-r^k F {0)-ir- 
' ' — \§(k = k F (9))\ 



(36) 



where k F {9) is the Fermi wave number of the polar co- 
ordinate 9 and d^/dk is the derivative with respect to 
k = |k|. If we expand z(9) as 



z(9) = Y. z " em9 > 

n 

the gap equation is given as 

n 

where Xmn are the matrix elements of x(#, 9'): 



2tt 



d9d9' PF {9')e- tme x(9,9')e 



I\„in8' 



(37) 



(38) 



(39) 



The number of basis functions kept in solving the 
eigenequation is 30 to 40 in this paper. The k-space 
is divided into 200 x 20 points on equally spaced mesh in 
the numerical calculations of x(9, 9'). 



FIG. 1: Square lattice with next-nearest transfer t' . 
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FIG. 2: The exponent X £LS £L function of the electron density 
for t' = . (See [3. We have included x for the E represen- 
tation.) Since the line for Ai mostly coincides with that for 
B2, the Ai line is omitted. 



B. Simple square lattice 

Let us investigate the phase diagram for the square 
lattice (Fig|TVX|. The basis functions {e mf> } (n = 
0, ±1, ±2, . . . ) are classified into irreducible representa- 
tions according to the symmetry group. The eigenfunc- 
tions are specified by one of irreducible representations 
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FIG. 3: The exponent i as a function of the electron density FIG. 5: The exponent a; as a function of the electron density 
fort' = -0.1. fort' = -0.3. 
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0.12 - 



0.08 



0.04 - 




FIG. 4: The exponent X ELS £L function of the electron density 
for t' = -0.2. 



FIG. 6: The exponent i of Bi symmetry as a function of the 
electron density for t' = 0, -0.1, -0.2, -0.3 and -0.4. 



of the square lattice (see Table 1). It is convenient to use 
real basis functions cos(n0) and sin(n0) for this purpose. 
The gap function in each representation is(l9j 



5(0) = z«cos(4£6>) A x , 



e=x 



z{6) = ^ z«sin(4^0) A 2 , 



(40) 



(41) 



r(0) = ^z«_ 2 cos(4£- 2)0 Bx, (42) 



r(0) = ^ z«_ 2 sin(4^ - 2)0 S 2 . (43) 



z(0) = ^ z 2 £-iCos(2^ - 1)0 E. 
e=i 

or sin(2^ - 1)0 (44) 

In Ref.[19j] the representations Ai~B2 were investigated. 
Here the E symmetry for triplet pairing is also examined. 
The eigenequation is solved for the above shown basis 
functions in the space of each irreducible representation. 
The eigenvalue x for t' = is shown in Fig UV Al as a 
function of the electron density n e . For n e > 0.6 the 
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FIG. 7: Phase diagram in the n e -t' plane for t' < 0. s 
denotes the pairing state with extended-s wave symmetry. In 
the s-wave region for small the s- and p-wave states are 
sometimes nearly degenerate. Small regions near boundaries 
are not shown. 
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t' 



FIG. 8: Phase diagram in the n e -t' plane for t' > 0. s, g 
and d x 2_ y 2 wave pairing states are almost degenerate in the 
low-carrier region for large t' . 



paired state with d x 2_ y 2 symmetry is most stable for t 1 — 
0. Since the exponent x sensitively depends on the van 
Hove singularity, x is an increasing function of n e near 
half filling for t' = 0. 

The exponent x for t' = —0.1, —0.2 and —0.3 is shown 
in Figs IIV A[ IIV Al and IIV Ai respectively. The expo- 
nents for small electron filling are not shown here because 
the high numerical accuracy is required for exponentially 
small exponents. As is shown in the figures, the d-wave 
state is most stable near half-filled case for t' in the range 





t 








A' 























FIG. 9: Square lattice with anisotropic next-nearest transfer 
t' (anisotropic triangular lattice) which is the lattice of organic 
conductors. 



x 




0.4 0.6 0.8 1 



n 

FIG. 10: The exponent x on the square lattice with 
anisotropic t' = —0.1. 



of < t' < 0.4. The position of the van Hove singular- 
ity depends on t', and the peak of x shifts as — t' > 
increases (Fig llV A| . x has a sharp peak showing a log- 
arithmic increase due to the van Hove singularity: 

x ~ -log|^i - fj, vH \, (45) 

where y, v n is the chemical potential corresponding to 
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FIG. 11: The exponent x on the square lattice with 
anisotropic t' = —0.5. 
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FIG. 13: Phase diagram for the square lattice with 
anisotropic t' > (lattice of organic conductors). 
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however, indicated that the optimum doping rate is al- 
most the same for Bi2212 and LSCO^. This may be 
a flaw of the weak coupling formulation, which, however, 
may not be completely remedied by the strong coupling 
treatment since the van Hove singularity still has a large 
effect on the critical temperature. This suggests that we 
must reexamine the structure of the Fermi surface of high 
temperature cuprates. In particular, the band parame- 
ters for Bi2212 will be modified if we take into account 
the double layer structure 27, 28]. The band structure re- 
ported by recent studies U 2°| is well fitted using smaller 
t' such as (si 



t' ~ -0.2. 



(46) 



FIG. 12: Phase diagram for the square lattice with 
anisotropic t' < (lattice of organic conductors). 



the van Hove singularity. The figure suggests higher T c 
for small — t'. The antiferromagnetism, however, may 
compete and suppress superconductivity near half filling. 
Hence we must have a bell-shape critical temperature as 
a function of the electron filling. 

It was pointed out from the electronic states calcu- 
lations that the Fermi surface is much deformed for 
Tl 2 Ba 2 Cu0 6 ,[2| and HgBa 2 Cu0 4 [Il] for which the 
band parameter values must be assigned as t' ~ —0.4 and 
t" ~ 0.1 (third-neighbor transfer). Bi 2 Sr 2 CaCu 2 08+,5 
(Bi2212) also has deformed Fermi surface so that t' ~ 
—0.3 and t" ~ 0.2[25]. For these values the optimum 
doping rate must be larger than that for Lai-zSr^CuC^ 
(LSCO) for which t> 0.1 and t" ~ 0. Experiments, 



The phase diagram in the n e -t' plane is shown in 
FigUvS for t' < and in Fig. ITVAl for t' > 0. For 
n e ~ 0.5 and — t' ~ 0.4, there is a possibility that 
the p-wave superconductivity is realized. For ex amp le, 
the ruthenate superconducting material Sr 2 Ru04|31| is 
sometimes modeled by the one-band Hubbard model for 
the 7 orbital with t' ~ —0.4 and n e ~ 0.67 after the 
electron-hole transformation. The state of these param- 
eters just corresponds to the point within the singlet re- 
gion near the boundary to p-wave regions in Fig. IIV Al In 
order to obtain the stable p-wave pairing for the parame- 
ters corresponding to Sr 2 Ru04, we may need to consider 
the multi-band structure including a and (3 orbitals(32l|. 
For t' > we have a large d-wave region. 

If t' is large and negative, i.e. if — t' > 0.5, we have the 
case with two Fermi surfaces; one is a large Fermi sur- 
face (FS1) and the other is a small Fermi surface (FS2) 
inside of the larger one. In this case we must examine 
the coupled equation of two gap functions z\ and z\ cor- 
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TABLE I: Irreducible representations of C± v for the square 
lattice. One of basis functions are also shown. 



Rep. Symmetry Bases 



M 


s 1 




cos(46») 


A 2 


9 


xy(x 2 - y 2 ) 


sin(40) 


Bi 


d x 2_ y 2 cos{k x )-cos(k y ) 


cos(26>) 


B 2 


d xy sa\(k x )sm(k y ) 


xy 


sin(20) 


E 


p sm(k x ),sm(k y ) 


x,y 


cos(6»),sin(6») 



responding to two Fermi surfaces: 

A y, x 11 (k+k')4^(ao+| F E x 12 (k+k')4 = 

k'-.FSl k':F52 

(47) 

k':_FSl k':FS2 

(48) 

where the symbol X)fe' _FSi indicates the summation over 
the Fermi surface FSi and x y (k+k') is the susceptibility 
for k on FSi and k' on FSj. The stable pairing symmetry 
is also obtained using the electron-hole transformation for 
t' > for which we have almost only one Fermi surface 
even in the electron-doped case. 

C. Square lattice with anisotropic t' ^ 

i 

The Hubbard model on the square lattice with 
anisotropic next-nearest-neighbor transfer t' (Fig UV A|) 
has been investigated intensively as a model for org anic 
conductors such as BDET-TTF(ET) moleculesfll H 
35] The model for organic conductors is well known as 
the Hubbard model with anisotropic next-nearest neigh- 
bor transfer t' (which is sometimes called the anisotropic 
triangular lattice). The dispersion relation is 

£fc = — 2t(cosk x + cosky) — 2t'cos(k x + k y ) — [i, (49) 

This model has the two-fold rotational symmetry and we 
classify the irreducible representation using the C^v point 
group (Table 2). The exponent x is in Figs II V Al and II V A1 
as a function of the electron density n e for t' = — 0.1 and 
t' = —0.5, respectively. As apparent from the figures, 
the d-wave state is stable over the whole region, which 
is consistent with the FLEX prediction (36[. The phase 
diagram in the n e -t' plane is presented in Fig UV Al for 
t' < and in Fig UV Al for t' > 0. For this model we 
conclude that the d-wave pairing is stable over the whole 
range of parameters. 

] 
t 

D. Three-band d-p model ' 

i 

The formulation is also applied to the three-band 
model for the CuC>2 plane 37]. We are interested in 



TABLE II: Irreducible representations of C-iv for the square 
lattice with anisotropic next-nearest-neighbor transfer. 



Representation Symmetry Bases 


Ax 




cos(26») 


A 2 


dxy 


sin(2#) 


Bi 


Px X 


cos(0) 


B 2 


Py y 


sin(0) 



.V 




-0.8 -0.4 0.4 0.8 



n 

FIG. 14: i as a function of the carrier density n e for the 
square lattice d-p model: n e > for hole doping and n e < 
for electron doping. 



0.4 - 



-0.4 - 




-0.8 I 1 1 1 1 

0.1 0.2 0.3 0.4 

t 

PP 

FIG. 15: Phase diagram for the three-band d-p model in 
the plane of the carrier number n and t pp in the range of 
< t PP < 0.4. We set e p — = 2 and tap = 1. n = 
indicates the half filling, and the positive and negative n are 
for hole doping and electron doping, respectively. 
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the relation between the single-band Hubbard model and 
the three-band d-p model. The pairing symmetry in the 
electron-doped cuprates is still controversial between the 
d-wave and s-wave order parameter [H, [H, |4(j. The 
Hamiltonian is 

H dp = e d ^24 a d i<7 +ep^2(p\ +s . /2(T Pi +i /2a 

ia ia 
+ Pi+y/2crPi+y /2<t) 

+ tdp '^^[d\ (J (pi+x/2a +Pi+y/2a —Pi-x/2cr —Pi-y/2a 
ia 

+ h.C.] + t pp ^2\p i+y -/ 2a Pi+x/2a - v\ + y/2ePi-&/2a 
ia 

~ Pl-y/2<rPi+£/2° + Pl-y/2aPi-i/^ + h - c -] 

+ u d J2 d h d n4i d n- ( 5 °) 

i 

In this subsection the energy is measured in units of t dp . 
The energy levels of the non-interacting Hamiltonian is 
written as in a concise form [53]: 

2 /0 k + 27raA e d -e p 

e - = vi £kcos (v — s — ; + (5l) 

for a = 0, 1 and 2, where 

*k = yjK\ 2 + \vZ\ 2 + (T&)* + (e d -e p )*/3, (52) 

4>k = | + sign{ Sii ) - arctaa v /|l-4te/(27a£)|^ , 

(53) 

su = (e d -e p ) - | + (<) 2 ) +<(«*+«)- 

(54) 

where ry^ = 2it dp sin(k x /2), rfc — 2it dp sin(k y /2), and 
77k = — 4tp J ,sin(fe x /2)sin(A^ / /2). e£ for a — 0,1,2 is the 
dispersion relation of the upper, lower and middle band, 
respectively. We examine the doped case within the hole 
picture where the lowest band is occupied up to the Fermi 
energy /z. The effective interaction is 

VW = ^ + f X dd (k + k'), (55) 

where 

p a/3 £ P £ q+P 

Here f£ is the Fermi distribution function, 

/-=(e« e £-") + l)~ 1 . (57) 
The weighting factor of d electrons is defined as 

(«£-«)(<£-«Z) 



where a, f3 and 7 are different from each other. The gap 
equation is 

A k = -J2™kV k d k d ,w k ,A k ,^-, (59) 



where w k = and E k = + A| for the lowest-band 
dispersion = ej. — /i. 

e?-wave pairing is predominant over the whole range in 
the parameter space as is shown in Fig UV Cl In partic- 
ular, d 2 ,2_ J/ 2-wave pairing is stable near half-filling. Al- 
though the extended s-wave pairing is possible in the nar- 
row region near half filling in the Gutzwiller variational 
Monte Carlo study [14|, we have no chance of s-wave su- 
perconductivity within the weak-coupling perturbation 
theory. The phase diagram for the d-p model is shown in 

Fignsia 



V. SUMMARY 

We have examined the phase diagram with respect to 
pairing symmetry on the basis of the two-dimensional 
Hubbard model. The weak coupling formulation is con- 
venient to investigate the phase diagram in detail. The 
results are almost consistent with the strong-coupling 
perturbation theory. We summarize the results as fol- 
lows. 

(1) The d-wave pairing is stable near half filling for the 
square lattice and the anisotropic square lattice. 

(2) The gap function has a maximum at the van Hove 
singularity. As the second neighbor transfer t' increases, 
the energy of the van Hove singularity decreases. For 
large t' = —0.3 ~ —0.4, the optimal doping is more than 
25 percent doping, i.e. n e < 0.75. For small third neigh- 
bor transfer t" the situation remains the same. The large 
—t' is assigned to several high-temperature cuprates to fit 
the angle resolved photoemssion spectroscopy (ARPES) 
data or the Fermi surface obtained by the band struc- 
ture calculations. Most of them, however, have optimum 
critical temperature in the range of 0.8 < n e < 0.85. 
Thus the weak coupling analysis suggests that we must 
consider other electronic or lattice interactions, or reex- 
amine the band parameters t' and t". Recent ARPES 
studies have reported the band structure which is well 
fitted using rather smaller t' such as t' ~ —0.2 by our 
analysis. 

(3) The predictions of the weak-coupling theory are 
almost consistent with the variational Monte Carlo 
method. An effective interaction to induce superconduc- 
tivity is possibly the simple x( c l) with renormalization in 
the Gutzwiller variational theory. 

(4) For the d-p model, the d-wave pairing is predominant 
in the wide range and the phase diagram is almost sym- 
metric between electron and hole dopings. Although the 
pairing symmetry in the electron-doped cuprates is con- 
troversial, only the <i-wave pairing is possible near half- 
filling in the weak-coupling formulation. 
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£&) 1 . Since we are interested in the third-order contri- 
butions, £„ (of the order of U 2 ) is neglected as follows: 



FIG. 16: The exponent x for the second-order (open symbols) 
and third-order (solid symbols) perturbation in U. We set 
U/t = 0.1. The symbol Non-V indicates the results obtained 
without vertex corrections. 
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APPENDIX A: HIGHER-ORDER CORRECTIONS 

In the Appendix we examine higher-order corrections 
to x. If the third-order terms have an effect to reduce 
the exponent x, the results obtained using the second- 
order perturbation have a possibility to become unstable 
as U increases. It is not an easy task to derive an effective 
Hamiltonian up to the third order of the interaction using 
the canonical transformation. The gap equation up to the 
third order of U has been obtained using the perturbative 
expansion for the Hubbard model[4ll |42|. The Green's 
functions satisfy the Dyson equations: 

G(k,ie„) = G (k, ie n ) + G (k, ie„)E„(k, ie„)G(k, ie n ) 
+ Go(k, ? )F*(k,ie n ), (Al) 



G(k,ie„ 



F{k, ien) 



4+^+|£ a (k,or 



(A3) 



el + e k + \X a (k,ien)r 1 ' 



The equation for the anomalous self-energy is 

£ a (k, ie n ) = — V [U + C/ 2 xo(k + k', ie n + ie n ,) 

+ 2U 3 Xo (k + k', ie n + ie n ,) 2 ]F{k', ie n >) 

+ U *Wn2 J2 Go(k',ie n ,)[xo(k + k',ie n + ie n ,) 

k' ,e n / ,p,e^ 

- 4> Q (k + k', ie n + ie„')]G (k + k' + p, ie n + ie n > + ie e ) 
x F(p,iet) 

+ U3 02pp J2 G'o(k',ie„0[xo(-k + k',-ie„ + i e ri 

k' ,e n / ,p,e£ 

- o (-k + k',-ie n + ie n >)] 

x G (-k + k'-p, — ie n +ie n > - iee)F(p,iei), (A5) 

for [3 = 1/kgT. The second and third terms originate 
from the vertex corrections. Xo(q.)^to) an d 0o(Q>^m) 
are defined as 



Xo(q,«Wm) = — > —7 -r^- 

N uv m + k,k- t,k+q 



4>o (q, iw r , 



m) - n-z-k +g ) 

Zk — £-k+q 



(A6) 



(A7) 



where u m = 2i:mkBT. We assume that E a is small and 
that we can neglect the e-dependence since we consider 
the small-?/ limit. We set A^ = E a (k, e n = 0), then 
the equation for Ak is derived. We show the results in 
FigfA]for U/t = 0.1 on the square lattice. The exponent 
x slightly decreases due to the third-order corrections. 
There is an cancellation among the third order terms. 
As has been shown in the literature 4l| , the vertex 
corrections reduce the exponent x and T c compared to 
that without vertex corrections. 



F(k,ie n ) = G (k (k, ie n )F(k, ie n ) 

- G (k, ze„)E a (k, ie n )G{-k, -ie n ), (A2) 

where e n = (2n + l^ksT is the Matsubara frequency, 
and £ n (E Q ) is the normal (anomalous) self-energy. Go 
is the free-electron Green's function: Go(k, ie n ) = (*£« — 
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